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Given a set of collective variables, a method is proposed to obtain the associated conjugated 
collective momenta and masses starting from a microscopic time-dependent mean-field theory. The 
construction of pairs of conjugated variables is the first step to bridge microscopic and macroscopic 
approaches. The method is versatile and can be applied to study a large class of nuclear pro¬ 
cesses. An illustration is given here with the fission of 258 Fm. Using the quadrupole moment and 
eventually higher-order multipole moments, the associated collective masses are estimated along the 
microscopic mean-field evolution. When more than one collective variable are considered, it is shown 
that the off-diagonal matrix elements of the inertia play a crucial role. Using the information on 
the quadrupole moment and associated momentum, the collective evolution is studied. It is shown 
that dynamical effects beyond the adiabatic limit are important. Nuclei formed after fission tend to 
stick together for longer time leading to a dynamical scission point at larger distance between nuclei 
compared to the one anticipated from the adiabatic energy landscape. The effective nucleus-nucleus 
potential felt by the emitted nuclei is finally extracted. 

PACS numbers: 24.75.-K, 21.60.Jz ,27.90.+b 


I. INTRODUCTION 

Nuclear time-dependent mean-field based on the en¬ 
ergy density functional approach is experiencing nowa¬ 
days a renewal of interest |lH6|. In particular, it allows 
one to describe a wide variety of dynamical processes 
ranging from small to large amplitude collective motions, 
including nuclear reactions. Among the most difficult 
challenges, one can mention the description of nuclear 
fission where a single nucleus encounters large deforma¬ 
tion leading finally to separated fragments. 

Although the fission process is quite well understood 
phenomenologically (tJ, it still remains one of the most 
difficult process to describe microscopically. One of the 
major difficulty is the necessity to treat collective and 
single-particle degrees of freedom (DOF) simultaneously 
as quantum objects @j. Moreover, fission is a dynami¬ 
cal process and therefore should be treated as such. To 
describe the dynamic of fission, one can a priori use two 
strategies. Using the fact that the time-scale associated 
to fission is rather large, the most common starting point 
is to first select a few collective DOFs and generate an 
adiabatic energy landscape. Then, the time-scale associ¬ 
ated to fission can be evaluated using semi-classical ap¬ 
proximation. Alternatively, the real time dynamic can be 
explicitly followed using for instance the Time-dependent 
Generator coordinates method |(|i|. Two important prob¬ 
lems are generally encountered in this strategy. First, 
this approach can hardly accommodate too many col¬ 
lective DOFs due to the increasing complexity. Second, 
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while rather slow, the dynamic of fission might deviate 
from the completely adiabatic path when the two frag¬ 
ments approach the scission. 

An alternative strategy is to use microscopic quan¬ 
tum transport theories such as the time-dependent en¬ 
ergy density functional (TD-EDF) theory. This approach 
offers the possibility to describe some aspects of the fis¬ 
sion process without assuming adiabaticity while leaving 
the possibility to explore rather complex shapes during 
the separation process. In addition, the recent inclusion 
of dynamical pairing has opened new perspectives E3- 
Il5j |. On the other hand, the TD-EDF approach cannot 
describe completely the fission process due to the absence 
of spontaneous symmetry breaking and due to poor treat¬ 
ment of quantal effects in collective space. However, it 
can still provide important information in particular af¬ 
ter the system has passed the fission barrier. This has 
been illustrated in Refs. iTol HU and more recently in 
Ref. 0. 

The TD-EDF directly performs the evolution of single¬ 
particle states in a self-consistent mean-field. From this 
evolution, one can directly infer the information on any 
one-body degree of freedom like multipole deformation, 
neck formation and/or fragment separation, final kinetic 
energies ... The aim of the present work is to explore the 
possibility to get macroscopic transport coefficients, like 
collective mass, collective potential or energy dissipation, 
directly from TD-EDF. In particular, this should allow 
us to compare the result of TD-EDF to similar quantities 
generally obtained in the adiabatic limit and/or macro¬ 
scopic models. 

In the following, we analyze first how collective masses 
and momenta can be associated to given collective ob¬ 
servables along a microscopic mean-field path. Once the 
pairs of conjugated operators are available, a macroscopic 
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reduction of the microscopic approach can be made to 
give physical insights. For the fission process, the en¬ 
ergy sharing between internal and selected DOFs can be 
precisely scrutinized. 


II. COLLECTIVE MASS AND MOMENTUM 
EXTRACTED FROM DYNAMICAL 
MEAN-FIELD THEORY 


In the present section, we assume that the mean-field 
trajectory including or not pairing is known, leading to 
a specific trajectory in the Liouville space of the nor¬ 
mal and anomalous density matrices (p(t), n(t)). Starting 
from these densities, we want to extract information on a 
set of given collective variables. Note that, although we 
present examples specifically on the fission process, the 
approach developed here is general and can be applied to 
other processes. 

Let us consider a given collective DOF associated to 
the one-body operator Q a . We restrict the present dis¬ 
cussion to the case where Q a corresponds to a local one- 
body operator, i.e. 

Qa = J d 3 r Q a (r)^(r)^(r), (1) 


where 'L'l(r) and 'I'(t’) are the creation/annihilation op¬ 
erators of a particle at position r. For the sake of simplic¬ 
ity, we omitted the spin and isospin quantum numbers. 
Note that most macroscopic DOFs of interest like mul¬ 
tipole operators, relative distance, mass asymmetry, ... 
correspond to expectation values of local one-body oper¬ 
ators. 

Along the mean-field trajectory, the collective evolu¬ 
tion is given by: 


q a (t) = Tr (Q a p(t)) 


Qa(r)n(r, t)d 3 r, (2) 


where we have introduced the local density n(r,f) = 
( r\p(t)\r). The first step to bridge the microscopic mean- 
field theory and a macroscopic like evolution for the col¬ 
lective variable q a is to find the corresponding conjugated 
momentum p a and associated collective mass, denoted by 
M a . Assuming a classical equation of motion for q a , one 
should fulfill the simple constraint q a = p a /M a . In the 
TD-EDF approach, we have: 

ih d (Qa) = TV (Q a [h(p),p \), 
dt 

where h{p) is the time-dependent mean-field hamilto- 
nian. Using properties of the trace and assuming a local 
mean-field potential, we obtain: 


d(Q a ) 

dt 


i 

2 hm 


Tr([<3 Q ,p 2 ] p(t)) , 


(3) 


where m is the nucleon mass. From this expression, a 
simple guess for the conjugated momentum is to assume 


that it can directly be defined through: 

Pa = II(*![£«,P 2 ]|j)a|aj> (4) 

ij 


where (aj,Gq) are creation/operators of any complete 
single-particle basis. From this definition, one can easily 
see that the matrix elements of the operator P a are given 
by: 


P n = -ih- 


M a f V 2 Q 0 


m 


VQ q .V 


(5) 


that is similar to the expression obtained in [20] using a 
variational principle around a static mean-field to study 
anharmonic effects in giant resonances. 

One shortcoming of the above expression is that the 
operator P a contains the collective mass M a that is un¬ 
known. To further progress, we seek for the condition 
that ( Q a ,Pa ) are conjugated observables. Similarly to 
the Time-Dependent RPA (TDRPA), we need to impose 
the condition: 


Tr (p(f)[Q a , P a ]J = ih , 


( 6 ) 


along the trajectory. We have: 

<[&»,£*]) = Tr ([Q a ,P a ]p{t)) = +ih—(VQ a .VQ a ) 

m 

We see that the condition (0 determines uniquely the 
collective mass through: 

i Tj = U WW ..V«.] 1 (7) 


all along the trajectory and henceforth also leads to an 
unambiguous definition of the collective momentum P a 
when reporting the mass in Eq. ©• Similar formula is 
sometimes used to compute collective mass from a micro¬ 
scopic adiabatic energy landscape (see for instance gH- 
[23jl. The difference is that this expression has been de¬ 
rived here without assuming adiabaticity. In addition, 
since the expectation value is directly performed using 
the time-dependent mean-field density it automatically 
contains possible influence of other DOFs as well as the 
pairing effects. 

The mass formula © is rather straightforward to cal¬ 
culate. We give illustration of some expression obtained 
for specific collective operators in Appendix [A] 

Once the mass and the momentum are known, one can 
also define the collective kinetic energy corresponding to 
the selected variable as 

*?,„«) = H = (8) 


A. Generalization to several collective degrees of 
freedom. 

Let us now consider a more general case where a set 
of N collective DOFs {Qa]a=\,N are selected. A naive 
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generalization to previous section is to assign to each vari¬ 
able Q a , a collective momentum P a with matrix elements 
given by Eq. ©. One should a priori also generalize the 
commutation relation ©. Using Eq. ©, we have: 

([Q a ,Pp})=ih^-(S7Q a -VQp) = ih^-, (9) 

m M a p 

where the off-diagonal mass matrix elements reads 

—i-- = -Tr [pWQc.VQf,]. (10) 

M a p[t) m 

This expression naturally extend the previous case and 
was also given in Ref. [ 23 . As shown in Appendix [B] the 
diagonalization of the mass gives new canonical pairs of 
collective operators ( Q' k ,Pk )> whose commutation rules 
identifies with the TDRPA ones and are given by 

([Q' k ,Pl]) = ihS kl . (11) 


The diagonalization is equivalent also to removing the 
correlation among the variables like multipole moments. 

These new operators are particularly useful to get sim¬ 
ple expressions for the evolution and collective energy. In 
particular, we have Q' a = P ! a /M’ a while the collective ki¬ 
netic energy is simply given by 


E {a} = P’lit) 

kin Z^2M'(f)' 


( 12 ) 


Once the set of collective variables is properly de¬ 
fined, macroscopic analysis of TD-EDF evolution can be 
made. Such a connection from the microscopic level to 
the macroscopic one is illustrated below for the fission 
process. 


III. APPLICATION TO THE FISSION OF 258 FM 

To illustrate the method presented in the previous sec¬ 
tion, we consider the case of 258 Fm that was the sub¬ 
ject of the recent work [Tol j . This nucleus is anticipated 
to have three different paths towards fission. In this 
work, we concentrate on the so-called symmetric com¬ 
pact shape. The energy landscape is obtained using the 
EV8 program with a constraint of the quadrupole mo¬ 
ment J24J. We use here the standard definition for mul¬ 
tipole moments 


^ = \[wfp* Y ^ < 13 > 

leading for instance to Q 2 {f) = 2z 2 — x 2 — y 2 . 

An illustration of the potential energy curve (PEC) is 
shown in Fig. [T] as a function of the quadrupole moment 
Q 2 . As in ref. [ij|, the Sly4d Skyrme functional [lj is 
used for the mean-field channel while a constant interac¬ 
tion is retained for the pairing channel. The static calcu¬ 
lations are performed with a mesh size 13.2 x 24.4 x 13.2 
fm 3 and a mesh step Ax = 0.8 fm. 
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FIG. 1. (Color online) Potential energy curve of 258 Fm nu¬ 
cleus as a function of the quadrupole deformation parameter 
(in barn unit). Isosurfaces of the total density drawn at half 
the maximum value at Q 2 = 34 b, 80 b, 194 b, and 399 b are 
also shown. The horizontal lines indicate the different start¬ 
ing points that are used in this work as initial conditions for 
the time-dependent evolution. The different vertical dashed 
lines corresponds from left to right to Q 2 = 34.2 b, Q 2 = 80 b 
(barrier position), Q 2 = 160 b (spontaneous fission threshold 
Qt) , O 2 = 182 b, Q 2 = 194 b, Q 2 = 296 b, and Q 2 = 400 b. 
The two thick arrows indicate the spontaneous fission thresh¬ 
old Q 2 h and the adiabatic scission point Q| c . 


The dynamical evolution of the system starting from 
any point of the PEC can be made consistently using 
the recently developed TD-EDF code including pair¬ 
ing in the BCS approximation [TEI . I 25 I [2f|. Dynami¬ 
cal calculations shown here are performed in a mesh of 
size 26.4 x 72.8 x 13.2 fm 3 with the same mesh step 
as in the static case Ax = 0.8 fm. The time step is 
At = 1.5 x 10 -24 sec ~ 0.45 fm/c. In the present cal¬ 
culations, reflection and axial symmetries are assumed 
in the constrained calculation. Since symmetry cannot 
be broken spontaneously by mean-field, only even multi¬ 
pole moments can be non-zero during the evolution. In 
particular, we do not consider here possible octupole de¬ 
formation. 

As it was observed previously including or not pair¬ 
ing, the system will spontaneously separate into two frag¬ 
ments only above a certain value of the initial quadrupole 
moment, which is larger than that of the fission barrier 
shown in Fig. [T] jl6Hl9l ]. The lowest initial quadrupole 
moment leading to spontaneous fission within TD-EDF is 
called hereafter “Dynamical Fission Threshold” and will 
be denoted by Q^ 1 . In the present calculation, the thresh¬ 
old deformation is approximately — 160 b. The 
shaded area in this figure indicates the region where the 
system does not spontaneously fission. The fact that Q 2 
is well beyond the expected barrier position signs the de¬ 
viation from the adiabatic limit of the microscopic trans¬ 
port theory close to single-particle levels crossing. This 
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FIG. 2. (Color online) Single-particle energies in 258 Fm nu¬ 
cleus along the adiabatic PEC. The green solid curves show 
the neutron and proton Fermi energies. Positive and negative 
parity states are respectively shown with red filled and blue 
open triangles. The vertical lines indicate the initial values of 
the quadrupole moment taken in the present dynamical cal¬ 
culations. The shaded area presenting the region where the 
system does not spontaneously fission is also shown. 


point was already discussed in Ref. [16). To illustrate 
the connection between the dynamical fission threshold 
and level crossing, the single-particle energies evolution 
obtained in the static constrained mean-field are shown 
in Fig. [2] as a function of Q 2 . 

We see in particular that for large Q 2 a gap in single¬ 
particle energies appears that signs the transition from 
one to two nuclei. At low quadrupole moments, many 
crossing occurs. When one leaves the system initially in 
the shaded area, single-particle wave functions will evolve 
in time. However, the motion is not adiabatic and both 
occupations above and below the Fermi energy will be 
populated in time. The PEC is meaningful only if low¬ 
est levels are preferentially occupied during the evolution 
while higher levels are depopulated. TD-EDF including 
or not pairing does not lead to sufficient reorganization 
of single-particle occupation numbers during the cross¬ 
ing to follow the adiabatic PEC as it has already been 
realized in Refs. flfH - Hfij . As studied in Ref. [l3], the 
system initialized inside the shaded area can still fission 
if, for instance, a boost in the quadrupole momentum is 
applied at initial time. In practice, the boost is imposed 
by applying the local operator exp(ip 2 Q 2 (r)/h) to each 
single-particle wave-function. This induces an additional 
initial collective kinetic energy E3 = 

ET = ^J |VQ 2 (r)| 2 n(r, t = 0 )d 3 r 

where n(r,t = 0) denotes the local density of the system 
in the adiabatic curve selected at a given initial moment. 


Note that E= p\/2M 2 with the quadrupole mass given 
by Eq. Q). 

In Fig. [3] several examples of density evolution ob¬ 
tained for different initial conditions including or not a 
boost initially and leading to fission are illustrated. We 
see that a variety of phenomena including ternary fission 
in some cases can be observed. 


A. Mass parameter from TD-EDF 


In the present section, we consider different initial 
quadrupole deformations between the fission barrier and 
the scission point. The scission point corresponding to a 
quadrupole deformation Q | c can already be seen in the 
shown in Fig. [Q It corresponds to the kink in the PEC 
appearing at Q 2 — 270 b. After the scission point, the 
PEC is nearly dominated by the Coulomb interaction be¬ 
tween the two fragments (See also Fig. [5]). 

As an illustration, we consider that the initial state cor¬ 
responds to Q Jf 1 = 160 b, that is a situation just above 
the spontaneous fission threshold. This initial condition 
is similar to the one considered in Ref. [l9). In particular, 
it has been shown that if the system is left initially with 
zero collective energy, the total final kinetic energy of 
fragments after TD-EDF evolution is compatible with ex¬ 
perimental observation. To study possible non-adiabatic 
effect, initial conditions with boost of varying intensity 
(including no boost at all) are used. 

The quadrupole moment is selected as the most rele¬ 
vant collective DOF. Using Eq. Oh the associated col¬ 
lective mass is then given by (appendix [Aj 

A = i (2(^(0) + <%«))) • (14) 

where (r 2 (f)) and ( Q 2 (t )) are respectively the root-mean- 
square radius and quadrupole moment along the path. 

To get physical insight it is interesting to consider the 
situation where the system is already about to get sep¬ 
arated into two fragments with a neck. Assuming sim¬ 
ply that the neck position is at the center of the whole 
system, quantities like mass, position, momentum, and 
intrinsic deformations of each fragment can be estimated 
through: 


(X(f))[i] = J X(r)n(r)Q(z)d 3 r , 
(*(*)>[ 2 ] = J X(r)n(r) [1 - 0(z)] d 3 r , 


where X is the local operator corresponding to the spe¬ 
cific quantity under interest, 0(z) is the Heaviside step 
function and [i = 1,2] is a label of fragments. For a 
di-nuclear system, the quadrupole mass can be recast as: 


1 

M2 


M*) 

m 2 


i?2 W + ^ E [2(r 2 (t)) [i] + (Q 2 (t)) W ] 

i= 1,2 
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(a) Q 2 ' = 34.2 b 
E 2 ' = 891.8 MeV 
t = 0.00 fm/c 



t= 143.90 fm/c 




t = 287.80 fm/c 



t = 359.75 fm/c 


(b) Q 2 ' = 79.8 b 
E 2 ' = 605.2 MeV 
t = 0.00 fm/c 


(c)Qj‘= 159.6 b 
E 2 ' = 0.0 MeV 
t = 0.00 fm/c 





t = 323.78 fm/c 


t = 485.66 fm/c 



t = 431.70 fm/c 


t = 647.55 fm/c 



t = 539.63 fm/c 


t = 809.44 fm/c 



t = 431.70 fm/c 



t = 647.55 fm/c 


t = 971.33 fm/c 



t = 503.65 fm/c 



t = 755.48 fm/c 



t= 1133.22 fm/c 



(d) Qg' = 193.7 b 
E 2 ' = 0.0 MeV 
t = 0.00 fm/c 



t = 107.93 fm/c 



t = 215.85 fm/c 



t = 323.78 fm/c 



t = 431.70 fm/c 



t = 539.63 fm/c 



t = 647.55 fm/c 



t = 755.48 fm/c 



FIG. 3. (Color online) Density profiles obtained for different initial Qi. The densities are shown as a function of time. For 
initial Q 2 below the fission threshold, a quadrupole boost has been imposed initially. From left to right, the initial Q 2 values 
are Q 2 = 34.2 b, Q 2 = 80 b, Q 2 = 160 b and Q 2 = 194 b. The system is eventually initially boosted leading to non zero values 
of E 2 11 directly indicated in the figure. The iso-density curves are drawn from 0.03 to 0.15 fm“ 3 with increment of 0.03 fm~ 3 . 
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FIG. 4. (Color online) Top: Quadrupole mass parameter 
calculated from TD-EDF paths. In all cases, the initial 
quadrupole moment is Q 2 " = 160 b. Different trajectories 
correspond to different initial boosts. The corresponding ini¬ 
tial collective energies E % 11 are systematically reported in the 
figure. The quadrupole mass obtained using Eq. 0 assuming 
that the system follow the adiabatic PEC is also shown (solid 
line) for comparison. Bottom: Quadrupole mass obtained 
with varying initial quadrupole deformation, Q 2 nl = 80 b, 114 
b, 160 bj and 194 b. 


where /x(t) = mA^A^jA is the reduce mass of the system 
and R(t) is the relative distance between the center of 
mass of the two fragments. 

At very large distance, R(t) —> + 00 , we see that the 
mass is dominated by the first term and tends to infinity. 
For display purpose we consider, as a reference mass, the 
mass obtained assuming no intrinsic quadrupole deforma¬ 
tion (tj) 2 (t)}[i— 1 , 2 ] = 0 and using the simple prescription 

(r 2 (t))[j] = |Ai ^1.2fm 2 . The reference mass ob¬ 
tained in this way is denoted by M 2as . In the following, 
the quadrupole mass will always be shown with respect 
to this mass. 

In Fig. QJ-a, the ratio M 2 /M 2 ia s of quadrupole mass 
deduced with the present method is shown as a function 
of Q 2 along the TD-EDF trajectories for Q™ = 160 b. 
In this figure, the mass obtained using increasing initial 
boost in the quadrupole collective momentum are shown. 
To illustrate the departure from the adiabatic path, we 
also show the result obtained for a given Q 2 assuming 
that the local density identifies with the corresponding 
density directly obtained from the constrained mean-field 
calculation. In the following, the latter is referred to as 
“static mass”. 

We observe in Fig. 0a that the mass is in general 


rather close to the static mass, especially if E™ 1 = 0 MeV. 
In that case, the system first follows closely the adiabatic 
case and then some deviation is observed. The deviation 
occurs around the scission point. At this point, the slope 
of the PEC suddenly change to match with the Coulomb 
case that dominates at large distance. This increase of 
slope is expected to induce also a larger collective veloc¬ 
ity and therefore also induce a possible departure from 
the adiabatic limit. We see in this figure that the mass 
also depends on the initial collective velocity imposed to 
the system. The larger is the initial velocity, the more 
deviation from the static mass is observed. It is how¬ 
ever worth mentioning that the adiabatic/non-adiabatic 
behavior cannot easily be concluded solely from the dif¬ 
ference of mass as will be further illustrated below. 


B. Mass parameter for Q 1 ™ < Q 2 h 

As we have mentioned already, TD-EDF cannot spon¬ 
taneously lead to the separation of the system into two 
fragments below Q 2 h due to the improper treatment of 
level crossing. Still it is possible to induce a fission by 
imposing some collective velocity initially. As already 
noted in Ref. the collective energy that should be 

initially put in the system to observe fission is in gen¬ 
eral rather large. This aspect is further illustrated in 
Fig. 0 where we investigated systematically the minimal 
collective energy necessary to induce fission for selected 
initial quadrupole deformation. To obtain this curve, for 
each Q we systematically performed TD-EDF calcula¬ 
tion by increasing progressively the boost intensity. The 
error bars correspond respectively to the largest (resp. 
lowest) collective energy where fission is (resp. is not) 
observed. Note that for the two lowest Q'J U no binary 
fission but ternary fission is observed (see also Fig. 0. 

The collective energy are very high compared to the 
typical barrier height to fission. We would like to mention 
that this is clearly a pathology of TD-EDF at small initial 
deformation and beyond mean-field effects should clearly 
be included to obtain meaningful information from mi¬ 
croscopic transport models around the fission barrier. 
Still, to illustrate that the present method can apply in 
situation rather far from the adiabatic limit, we deduced 
the mass parameter for such initially highly excited sys¬ 
tems. A few examples are shown in panel (b) of Fig. 
0 In that case, important deviation are observed in the 
mass parameters compared to the static/adiabatic limit. 
Since large collective velocities are imposed initially, such 
deviations are not surprising. However, the difference 
can also stem from the fact that the initial boost can 
induce a motion that is not described by the simple one¬ 
dimensional energy landscape shown in Fig. 0 In par¬ 
ticular, preparing the system using constraint mean-field 
+ boost is a rather arbitrary choice that will induce spe¬ 
cific motion not only in the Q 2 collective space but also 
in a larger space of collective variables like the monopole 
(r 2 ), hexadecapole Q 4 , .... In particular, since the mass 
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FIG. 5. Minimum collective energy after a boost that should 
be initially deposited into the system to induce the fission of 
258 Fm for Q™ < Q 2 h - This energy is plotted as a function of 
the initial quadrupole deformation considered. For each Q 2 , 
the highest (resp. lowest) value of the initial collective energy 
where the scission is not (resp. is) observed is reported. 
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reported in Fig. [4] are compared for the same Q 2 , differ¬ 
ences observed between the static and dynamical masses 
stem from differences in the root mean-square radii that 
ultimately come from the differences in local densities. 
Clear differences are observed between densities shown 
in Fig. [3] and densities of the adiabatic PEC (Fig. [TJ. 
The differences in local densities can of course come from 
non-adiabatic effects but also from a more complex path 
in a multi-dimensional potential energy landscape that 
could not be simply reduce to the ID picture of Fig. [T] 


C. Total versus collective kinetic energy 


As we mentioned in Sec. El the present method al¬ 
lows us to access the set of conjugate momenta and the 
collective kinetic energy (CKE) as well as the masses as¬ 
sociated with the set of collective coordinates. The CKE 
of the set of collective variables {Q q } can be obtained 
following section III Al diagonalizing the mass matrix and 
using Eq. m- Note that the diagonal and off-diagonal 
matrix elements of the mass are given explicitly in ap¬ 
pendix [A] 

In Fig. [6l the CKE associated to the quadrupole 
and/or hexadecapole moment are displayed as a function 
of time during the fission process. We also compare these 
energies to the total kinetic energy computed through 


77'tot _ 

-^kin — 


2m 


d 3 r 


n(r , t) ’ 


(15) 


where j(r,t) is the single-particle current. Two different 
initial conditions are considered, one starting from an 
already elongated shape without boost and one with a 
more compact shape but where a boost in quadrupole 
momentum is applied to induce fission. 


FIG. 6. (Color online) Evolution of the total collective kinetic 
energy -Ekin as a function of time. The initial systems corre¬ 
spond respectively to a quadrupole moment (a) Q-? 1 = 171.0 
b without boost or (b) to Q 2 11 = 79.8 b with initial boost. In 
both cases, the CKE obtained using Eq. © and associated 
to Q 2 only (EjL), Q 4 only (Ej) in ) or both (E 2 ; + 4 ) are also 
shown. 


From this figure several interesting aspects could be 
seen: 

• At initial time Ej£* = E 2 in . This is indeed due 
to the fact that either the two are equal to zero 
(Fig. [HKa)) or that the initial condition (Fig. O-b) 
is such that all initial kinetic energy is imposed by 
the quadrupole boost. 

• The CKE associated to Q 4 is also initially non zero. 
This stems from the fact that Q 2 and Q 4 are not in¬ 
dependent collective variables. Therefore boosting 
in quadrupole moment also induces an excitation of 
the hexadecapole and most probably higher order 
even multipole moments. 

• Due to the rather strong correlations between Q 2 
and Q 4 , the off diagonal matrix elements of the 
inertia play an important role. Indeed, neglecting 
this contribution would give: 

^km-^kin + ^kin- ( 16 ) 

However, summing directly these two energies 
would exceed the total kinetic energy that is an 
upper bound whatever is the selected set of collec¬ 
tive variables. In Fig. [Gl accounts for the 

off-diagonal inertia and finally leads to an energy 
that is lower than E Jr*. 




















as a simple dissipative equation of motion: 



Q 2 (b) 


FIG. 7. (Color online) Evolution of the collective momentum 
as a function of time for different initial quadrupole defor¬ 
mations (with Q 2 m > Qf). The arrow indicates the scission 
point associated to the adiabatic potential. 


• At large distances, we see that 

ZT’tot ZT'2+4 7 t , 2 /-i 7 \ 

-^kin — -^kin ~ -^kin’ ) 

This is due to the fact that all kinetic energies are 
dominated by the relative motion of the two frag¬ 
ments in the exit channel. 


D. Collective evolution close to scission 

Here we investigate the collective evolution close to the 
scission point. The scission can be seen directly on Fig. 
[T] by the change of slope around Q| c ~ 270 b. The evo¬ 
lution of the collective momentum is displayed in Fig. [7] 
as a function of Q 2 for different initial deformations. We 
clearly see a different behavior depending if the initial 
quadrupole moment is above or below Q | c . As we will 
see below, for Q 1 ™ > Q s 2 the momentum evolution corre¬ 
sponds essentially to the evolution of two escaping nuclei 
boosted by their mutual Coulomb field. For Q If 1 < (5| c , 
the nuclear interaction between nuclei still plays a sig¬ 
nificant role and a richer evolution is seen. In that case, 
independently of the initial Q™ value, after some tran¬ 
sition time, all curves become nearly identical with one 
another. 

In the absence of dissipation and assuming that the 
dynamics stem uniquely from a collective potential, one 
would expect that the smaller is Q 1 ™, the higher is P 2 (t) 
as a function of <52- However, it clearly seems from Fig. 
[T]that part of the energy is dissipated in the early stage of 
the evolution. To further progress, we may follow Ref. @ 
and assume that the momentum evolution can be written 


P 2 = - 


dV co ii 

~dQ7 


1 dM 2 - 2 

2 8Q 2 


l(Q2)Q2, 


(18) 


where the collective potential V co u and the friction co¬ 
efficient 7 are a priori unknown quantities. In the adia¬ 
batic limit, the collective potential identifies with the one 
shown in Fig. [T|and y(Q 2 ) = 0 along the path. 

To remove the possible effect of the mass evolution and 
eventually access the potential and dissipative collective 
properties, it is convenient to define the quantity 

(19) 

This function is shown in Fig. [ 8 ] as a function of Q 2 for 
some of the evolutions presented in Fig. 0 If the macro¬ 
scopic transport equation m is valid, this quantity is 
expected to identify with: 

F(Q2) = -^ i -l(Q2)Q2. ( 20 ) 

0 Q 2 


and therefore is sensitive to both the potential and dis¬ 
sipative part. For comparison, we also display the cases 
where dissipation is assumed to be zero in Eq. (ESI) and 
where the potential part identifies either to the adiabatic 
potential or solely to the Coulomb field. In the latter 
case, the Coulomb potential for large relative distance, 
or large Q 2 , is approximated by 




ZiZ 2 e 2 

R 


Z^e 2 



1 

4 



The last expression is obtained for the symmetric fission 
case, i.e. Z\ = Z 2 = Z/2 and A\ = A 2 = A /2 consid¬ 
ered here, and further assuming no intrinsic quadrupole 
deformation of emitted fragments after scission. 

Fig. [8j-a gives interesting information on the different 
steps leading to fission. We first see that after the very 
first instant of the evolution where some dissipation oc¬ 
curred, all evolutions obtained with Q™ 1 fs Q| c are on 
top of each other. The dynamic before scission deviates 
significantly from the expected adiabatic one underlin¬ 
ing the importance of both non-adiabatic and dissipative 
effects. In particular, we clearly see that the dynami¬ 
cal formation of the neck differs from the adiabatic case. 
Defining the scission point as the Q 2 value where the neck 
density equal 1/10 of the saturation density, we observe 
that dynamically the scission occurs at much larger Q 2 
than the adiabatic case (arrows in Fig. 0. This has two 
consequences (i) first, the two nuclei stick together at 
larger distance compared to the adiabatic case. Accord¬ 
ingly, the nuclear field can play an enhanced role, (ii) we 
see that we should introduce the notion of “dynamical 
scission point” that a priori differs from the “adiabatic 
scission point” and that occurs at larger quadrupole mo¬ 
ment. In the present case of symmetric compact fission, 
the dynamical scission point occurs around Q s £ ,dyn ~ 350 
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FIG. 8. (Color online) (a) Function F(Q 2 ) obtained with TD- 
EDF using Eq. (1191) for the three evolutions with Q 1 ™ < Q 2 C 
displayed in Fig. [7] For comparison, we also show the forces 
acting on Q 2 that would be induced either by the adiabatic 
potential (black dashed line) or solely by the Coulomb field 
(black filled circle). The arrows in the figure indicate the 
Q 2 value where the neck density p nec k becomes ten times less 
than the saturation density p sa t = 0.16 fm -3 . The two arrows 
indicates the adiabatic path (static) and dynamical path (dy¬ 
namic). In the latter case, the position where p n eck/psat = 10 
is almost independent of Q'F- (b) Dynamical potential curve 
obtained by integrating F(Q 2 ) using Eq. (12311 . Again for 
comparison, the adiabatic potential and the Coulomb field 
are also shown. Here, the origin of energy is taken such that 
E = 0 for Q 2 —> 00 . 


b, compared to the adiabatic scission point Q 2 C,stat — 270 

b. 

After scission, the dynamical evolution is very close to 
the Coulomb field case (black filled circles). This indi¬ 
cates that no dissipation takes place after this point and 
that the dissipation mainly occurs at initial time of the 
calculation. For large Q 2 , we clearly observe some oscil¬ 
lations around the average Coulomb repulsion that could 
be attributed to the dynamical oscillation of the intrin¬ 
sic shapes of each nucleus. These oscillations obviously 
goes beyond the simple macroscopic approximation m 
since they involves additional intrinsic shape degrees of 
freedom. 

Following (gi], one could a priori use F(Q 2 ) to get the 
potential energy landscape as well as the friction coeffi¬ 
cient along the path. However, the method used in Ref. 
0 that consists in performing two evolutions with close 
initial conditions cannot be applied here due to the fact 


that the collective velocity becomes rapidly independent 
of Q™ 1 . Fig. [7] seems however to indicate that dissipa¬ 
tion occurs only at rather small Q 2 . For Q 2 > 300 b, 
one might assume that the motion is only driven by a 
potential denoted by V dyn (Q 2 ). Then, we have the ap¬ 
proximate relationship: 


dv dyn (Q 2 ) 

oq 2 


-F(Q 2 ) 


( 22 ) 


where F(Q 2 ) is estimated along the path using Eq. (fl9l) . 
To get the potential itself, one should fix the boundary 
condition. We know from Fig. [8j-a that the potential 
identifies to a good approximation with the Coulomb po¬ 
tential after the scission. The potential can eventually 
be obtained through the relation: 

gmax 

V dyn (Q 2 ) = V c (Qr x ) + [ 2 F{Q' 2 )dQ' 2 (23) 

JQ2 


where Q™ ax is taken much larger than the dynamical 
scission point. Note that we do not take into account 
here the excitation energy of the fragments. Examples 
of potential obtained in this way is shown in Fig. [Sj-b 
assuming Q™ ax = 800 b. 14a shown with the dashed 
curve in Fig. [8}b is drawn by shifting the PEC given 
in Fig. [T] so that it coincides with Vc at Q 2 = 433 b 
(> Q s 2 c ’ stat ) and thus I4d —> 0 for Q 2 —> 00 . 

We see in Fig. [SJ-b that the potential obtained using 
Eq. (1^1) differs significantly from the adiabatic one at 
small Q 2 due to dynamics and eventually non-adiabatic 
effects. Note that the dynamical potential should be in¬ 
terpreted with some caution since it might contain some 
dissipative effects especially at initial time. It is worth in 
particular mentioning that the adiabatic and dynamical 
potentials should be identical at initial Q 2 . We clearly 
observe in Fig. [BJb a lower value for the dynamical 
case. The difference between the adiabatic and dynami¬ 
cal curves at t = 0 corresponds to the energy transferred 
into the other collective degrees freedom or internal ex¬ 
citations during the fission. We see that this difference is 
« 23 MeV for Q'f 1 = 160 b. 


E. Dissipation estimated from energy balance 

Alternatively, in order to estimate the energy dissi¬ 
pated into the internal excitation of the fragments, an 
analysis similar to Ref. (id] has been made for the to¬ 
tal kinetic energy (TKE) of the outgoing fragments after 
fission. The TKE is defined as 

TKE = \^R 2 + - = T r ei + Vc{R) (24) 

Z H 

at large R well beyond the scission point. Note that the 
TKE is the energy of the relative motion between the 
fragments at infinite separation, and it is not the same 
quantity as either Eq. m or Eq. m , which may con¬ 
tain the internal excitation energy of the fragments. The 
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R (fm) 


FIG. 9. The relative kinetic energy T re i, the Coulomb energy 
Vc, and the total kinetic energy (TKE) of the fission frag¬ 
ments as a function of the relative distance R for Q 2 11 = 197 
b. 

energy dissipated into the other DOFs than the relative 
motion, or the excitation energy of fragments, is then 
given by E* = Eq — TKE, where Eq is the total energy 
of the system pjJ . Note that the origin of energy for Eq 
is taken as that for R(t) —> 00 with fragments staying at 
their ground states. 

In Fig. [9] we show relative kinetic energy T re[ , the 
Coulomb energy Vc, and the total kinetic energy (TKE) 
of the fission fragments as a function of the relative dis¬ 
tance R for the evolution with Q™ = 194 b. We see a 
plateau in TKE for large R at « 238 MeV, which is iden¬ 
tified as the TKE for this process, while the total energy 
of the system is Eq ss 256 MeV, which is given by the 
value of Kd in Fig. [BJb at Q™'- Taking the difference, we 
obtain E* « 18 MeV. This is close to the value obtained 
from the previous analysis confirming that the difference 
between the adiabatic potential and the dynamical po¬ 
tential estimated from Eq. (E3l) most likely stems from 
the energy dissipated along the fission path. 

IV. CONCLUSION 

In the present work, a method is proposed to construct 
conjugated collective momenta associated to a given set 
of local collective variables along a time-dependent EDF 
path. A detailed discussion is made on the proper defini¬ 
tion of associated inertia including the effect of its pos¬ 
sible off-diagonal matrix elements. Once pairs of conju¬ 
gated collective variables are obtained, one can make a 
macroscopic reduction of the microscopic mean-field dy¬ 


namic. 

An illustration is given here with the fission process. 
A precise analysis is made in the symmetric fission of 
258 Fm. The mass matrix is calculated along the fission 
path including only quadrupole moment, and/or both 
the quadrupole and hexadecapole moment. In particu¬ 
lar, the important role of off-diagonal matrix elements of 
the mass is underlined. Then, a detailed analysis of the 
macroscopic evolution in the quadrupole collective space 
is made. The importance of dissipation in the early stage 
of the evolution is discussed. Clear non-adiabatic effects 
are probed, in particular associated with the specific neck 
evolution. We show that the scission point, called dy¬ 
namical scission point, occurs at much larger quadrupole 
moment compared to that for the adiabatic path. An 
attempt to extract the nucleus-nucleus potential felt by 
the daughter nuclei after fission is also made. It is shown 
that this potential significantly differs from the adiabatic 
one due to the non-adiabatic effects and the dissipation of 
energy into the intrinsic excitations of fission fragments. 

The method presented here is rather versatile and 
could be used in other dynamical processes. For instance, 
it could a priori be used to study anharmonicity in giant 
resonances as well as possible coupling between collective 
modes. 
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Appendix A: General formula for the inertia matrix 

In this appendix we give the explicit expression for the 
inertia matrix 

7TL 

= Tr {pVQc-VQp], (Al) 

M a p 

associated with the general multipole moment 

&= \fl <a2) 

Its expressions for 2 < A < 6 are given in Table U| 

Using the Racah algebra technique, one obtains the 
following for the matrix 


m 

M\\t 


+ !)■(— 


(s - 2 L)!(s - 2A)!(s - 2A')!(s/2)!(s/2 - 1)! 

1)! (s/2 - L)!(s/2 -L- l)![(s/2 - A)!] 2 [(s/2 - A')!] 2 


•Tr [pr x+x '- 2 - L Q L \, 


(A3) 


where s = A + A' + L and the value of L runs over L = |A — A'|, |A — A'| + 2, |A — A'| + 4,..., A + A' — 2. It allows 
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TABLE I. Multipole operators Q\ = \J 2 \+i r> 'Y\o for A = 
2, 3, 4, 5,6. 


A 

operator 

2 

o 2 2 2 

2 z — x — y 

3 

2 z 3 — 3 zx 2 — 3y 2 z 

4 

i(35 z 4 — 30 r 2 z 2 + 3r 4 ) 

5 

i(63 z 5 — 70r 2 z 3 + 15 r 4 z) 

6 

i(231 z 6 - 315rV + 105rV - 5r 6 ) 

8 


TABLE II. Illustration of the inertia matrix elements m/M \ A / for the multipoles of A, A' = 2,3,4. The expressions for the 
multipole operators are given in Table [T] 


A \ A' 

2 

3 

4 

2 

4(2(r 2 ) + (Q 2 » 

y(3(r 2 Qi) + 2(Q 3 » 

|(9(r 2 Q 2 ) + 5(Q 4 » 

3 

— 

y(7(r 4 >+4(r 2 Q 2 } + 3(Q 4 )) 

Yf (99{r 4 Qi) + 77(r 2 Q 3 > + 150(Q 5 » 

4 

— 

— 

(462 (r 6 ) + 275(r 4 Q 2 > + 243<r 2 Q 4 ) + 175(Q 6 }) 


us to compute a matrix element for any multipolarities A and A'. In Table EH we illustrate expressions of the matrix 
elements for some important multipoles 2, 3, and 4, which are be related to elongation, mass asymmetry, and size of 
neck, respectively. 


Appendix B: Definition of new operators ( Q' a ,P ' a ) 


In this appendix, we give some intermediate steps to 
obtain operators that fulfill the commutation rules m 
along the TD-EDF trajectory. 

The momentum operator conjugated to the coordinate 
Q a is given by 


Pry = —ih 


Mry 


V 2 Q 0 


VQ a ■ V 


(Bl) 


the mass, i.e. 


a (3 ^ 


M nB Wf)l 6kl M' k 


(B3) 


where 1 /M' k are the eigenvalue of the inertia tensor. Then 
we have 


E WkaWgidQq, Pp/Mpp}) = — ■ (B4) 

a/3 


Introducing the new set of conjugated operators 


Accordingly, we have: 

([Q*,Pp]) =ih^-(S/Q a -VQp). 

m 

Introducing the mass matrix elements m/M a p = (’ VQ a ■ 
VQp): above expression can be written as 

([Qa>Pp/Mpp]) = iti-jrr —■ (B2) 


Q'k = E WkcQct, and P' =M'E W ka ^,(B5) 

a a. 

we see that these operators respect the desired commu¬ 
tation relation 

([Qlp;/Ml}}=ih6 kl ^r (B6) 

or equivalently 


We introduce the orthogonal matrix W that diagonalizes 


([Q'k: -P/]) — ifiSkl- 


(B7) 
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